Astronomy & Astrophysics manuscript no. romero 


February 1, 2008 


(DOI: will be inserted by hand later) 





o 
o 

>> 



6 



> 

oo 

OS 
(N 

vn 
o 
l> 
o 



X 



The formation of spiral arms and rings in barred galaxies 

M. Romero-Gomez^'^, E. Athanassoula\ JJ. Masdemont^ , and C. Garcia-Gomez^ 

* LAM, Observatoire Astronomique de Marseille Provence, CNRS, 2 Place Le Verrier, 13248 Marseille Cedex 04, France 
^ D.E.I.M., Universitat Rovira i Virgili, Campus Sescelades, Avd. dels Paisos Catalans 26, 43007 Tarragona, Spain 
^ I.E.E.C & Dep. Mat. Aplicada I, Universitat Politecnica de Catalunya, Diagonal 647, 08028 Barcelona, Spain 

Received 19 March 

Abstract. In this and in a previous paper (Romero-Gomez et al. 2006 i we propose a theory to explain the formation of both 
spirals and rings in barred galaxies using a common dynamical framework. It is based on the orbital motion driven by the 
unstable equilibrium points of the rotating bar potential. Thus, spirals, rings and pseudo-rings are related to the invariant 
manifolds associated to the periodic orbits around these equilibrium points. We examine the parameter space of three barred 
galaxy models and discuss the formation of the different morphological structures according to the properties of the bar model. 
We also study the influence of the shape of the rotation curve in the outer parts, by making families of models with rising, fiat, 
or falling rotation curves in the outer parts. The differences between spiral and ringed structures arise from differences in the 
dynamical parameters of the host galaxies. The results presented here will be discussed and compared with observations in a 
forthcoming paper. 
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1. Introduction 

Bars are a very common feature of disc galaxies. In a sam- 
ple of 186 spirals drawn from the Ohio State University Bright 
Spiral Galaxy Survey, Eskridge et al. ( 120001 ) find that in the 
near infrared 56% of the galaxies are strongly barred, while 
an additional 6% are weakly barred. Only 27% can be clas- 
sified as non-baiTed. A large fraction of barred galaxies show 
two clearly defined spiral arms (e.g. Elmegreen & Elmegreen 
1 19821 1. often departing from the end of the bar at nearly right 
angles. This is the case for instance in NGC 1300, NGC 1365 
and NGC 7552. Deep exposures, moreover, show that these 
arms wind around the bar structure and extend to large dis- 
tances from the centre (see for instance Sandage & Bedke 
I19941 I. Almost all researchers agree that spiral arms and rings 
are driven by the gravitational field of the galaxy (see Toomre 
119771 and Athanassoula 1984 , for reviews). In particular, spi- 
rals are believed to be density waves in a disc galaxy (Lindblad 
I1963I I. Toomre d 19691 1 found that the spiral waves propagate to- 
wards the principal Lindblad resonances of the galaxy, where 
they damp down, and thus concludes that long-lived spirals 
need some replenishment. 

There are essentially three different possibilities for a spiral 
wave to be replenished. First, it can be driven by a companion 
or satellite galaxy. A direct, relatively slow and close passage 
of another galaxy can form trailing shapes (e.g. Toomre [19691 
Toomre & Toomre [19721 Goldreich & Tremaine [T9781 [19791 
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Toomre [l 9 8 1 1 and references therein). They can also be excited 
by the presence of a bar. Several studies have shown that a ro- 
tating bar or oval can drive spirals (e.g. Lindblad 1960 ; Toomre 
[1969. Sanders & Huntley ^976^ Schwarz , 1979 . 1981; Huntley 
I1980I I. Athanassoula d 19801 1 studied the self-consistent quasi- 
stationary response of a galaxy composed of both gas and stars 
to a growing and rotating bar She showed that the response of 
both components is bar-like up to corotation, where it turns into 
a trailing two-armed spiral, ending approximately at the Outer 
Lindblad Resonance (OLR). The third alternative, proposed by 
Toomre ( |19811 l. is the swing amplification feedback cycle. This 
starts with a leading wave propagating from the center towards 
corotation. In doing so, it unwinds and then winds in the trail- 
ing sense, while being very strongly amplified. This trailing 
wave will propagate towards the center, while a further trail- 
ing wave is emitted at CR and propagates outwards, where it is 
dissipated at the OLR. The inwards propagating trailing wave, 
when reaching the center will reflect into a leading spiral which 
will propagate outwards towards the CR, thus closing the feed- 
back cycle. Note that, if there is an Inner Lindblad Resonance 
(ILR), the wave propagating inwards is damped at that reso- 
nance and the cycle is cut. 

Strongly baiTed galaxies can also show prominent and 
spectacular rings or partial rings. The origin of such morpholo- 
gies has been studied by Schwarz ( 119811 119841 [1985) . who 
followed the response of a gaseous disc galaxy to a bar per- 
turbation. He proposed that ring-like patterns are associated to 
the principal orbital resonances, namely ILR, CR, and OLR. 
The ILR would be responsible for the nuclear rings, CR would 
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be associated with the inner rings, which are indicated by an 
r in the de Vaucouleurs classification, and the OLR would be 
the origin of the outer rings, which are indicated by an R pre- 
ceding the Hubble type. Nuclear rings are small rings of star 
formation often found near the centres or nuclei of early-type 
barred galaxies, but we will not discuss these structures here. 
Inner rings are the well-defined rings that encircle the bars of 
barred galaxies. They are sometimes active sites of star forma- 
tion. Outer rings are the larger, more diffuse rings, about twice 
the size of the bar When these structures are incomplete, we 
designate them with the term pseudo-rings. There are different 
types of outer rings. Buta (1995) classified them according to 
the relative orientation of the ring and bar major axes. If these 
two axes are perpendicular, the shape of the ring is similar to 
an "8" and the outer ring is classified as /?i . If the two axes are 
parallel, the outer ring is classified as /?2- Finally, if both types 
of rings are present in the galaxy, the outer ring is classified as 
RiR2- 

In this paper, we develop further the basic idea that we pro- 
posed in a previous paper (Romero-Gomez et al. 120061 here- 
after Paper I). Rings and spiral arms are the result of the orbital 
motion driven by the invariant manifolds associated to periodic 
orbits around unstable equilibrium points. In Paper I we gave 
a complete and thorough description of the dynamical frame- 
work and, as an example, specifically calculated the orbital evo- 
lution of a particular type of rRi ringed galaxy. In this paper, 
we construct families of models based on simple, yet realistic, 
barred galaxy potentials. In each family, we vary one of the 
free parameters of the potential and keep the remaining fixed. 
For each model, we compute numerically the orbital structure 
associated to the invariant manifolds. In this way, we are able 
to study the influence of each model parameter on the global 
morphologies delineated by the invariant manifolds. 

In Sect.|2] we present the equations of motion and the galac- 
tic models that will be used for the computations. In Sect. [3] 
we give a brief description of the invariant manifolds and de- 
scribe their role in the transfer of matter within the galaxy. In 
Sect. |4] we show the different morphologies that result from the 
computations of each model and in Sect. |5]we briefly summa- 
rize. Applications to real galaxies, as well as comparisons with 
observations and with other theoretical work (Schwarz 1981 
[19841 [T9851 Kaufmann & Contopoulos ,1996 ; Patsis 2006; 
Voglis, Stavropoulos & Kalapotharakos l2006l l will be given in 
paper III of this series. 

2. Models 

2.1. Equations of motion and equilibrium points 

We model the potential of a barred galaxy as the superposition 
of three components, two of them axisymmetric and the third 
bar-like. This latter component rotates anti-clockwise with an- 

aular velocity Qp = QpZ, where Qp is a constant pattern speed 
. The equations of motion in this potential in a frame rotating 
with angular speed Hp in vector form are 

f = -VO - 2(np X f) - Qp X (Hp X r), (1) 

' Bold letters denote vector notation. The vector z is a unit vector. 



where the terms -2np x f and -Hp x (Hp x r) represent the 
Coriolis and the centrifugal forces, respectively, and r is the 
position vector. 

Following Binney & Tremaine (I1987I I, we take the dot 
product of Eq. ([T]i with f , and by rearranging the resulting equa- 
tion, we obtain 

dt 

where 

£,= i|r|2 + (i.-i|£ipxr|2. 

The quantity Ejisa constant of the motion and is known as the 
Jacobi integral, or as the Jacobi constant. Note that it is the sum 
of jT'^ + <t), which is the energy in a non-rotating frame, and of 
the quantity 

-^\£lpxrf^~£ll(x^+y\ 

which can be thought of as the "potential energy" to which the 
centrifugal "force" gives rise. Thus, if we define an efifective 
potential 

Eq. ([B becomes 

r = -V«>eff - 2(Qp X f), 
and the Jacobi constant is 

Ej^^\rf+ (Deff, 

which, being constant in time, can be considered as the energy 
in the rotating frame. 

Setting xi = x, X2 - y, xj, = x, and x^ - j, we can 
write the components of the equations of motion as 

' Xi = /i(xi,X2,X3,X4) = X3 

X2 = /2(xi,X2,X3,X4) = X4 

X-i — /3(xi, X2, X3, X4) = 2QpX4-<I)v 

X4 — /4(xi,X2,X3,X4) = — 2QpX3 — Oy, 

where we have defined O, = '^^^ and = Note that we 
restrict ourselves to the z - plane because the motion in the 
vertical direction essentially consists of an uncoupled harmonic 
oscillator and does not affect the motion in this plane (see Paper 
I). 

The surface Oetj - Ej is called the zero velocity surface, 
and its intersection with the z = plane gives the zero velocity 
curve. All regions in which <l)eft > Ej are forbidden to a star 
with this energy, and are thus called forbidden regions. The 
zero velocity curve also defines two different regions, namely, 
an exterior region and an interior one that contains the bar The 
interior and exterior regions are connected via the equilibrium 
points (see Fig. 2b of Paper I). 

The bar component is rotating as a solid body with a con- 
stant pattern speed Q.p. For our calculations we place ourselves 
in a frame of reference corotating with the bar, and the bar 
semi-major axis is located along the x axis. In this rotating 
frame we have five equilibrium Lagrangian points. Three of 
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these points are stable, namely L3, which is placed at the centre 
of the system, and L4 and L5, which are located symmetrically 
on the y axis. The two equilibrium points left, Li and L2, are 
unstable and are located symmetrically on the x axis. The po- 
sition of the corotation radius will be determined by the free 
parameter r/,, that is, the distance from the centre to Li, or, 
equivalently, to L2, because the model is symmetric with re- 
spect to both X and y axes. In our case, we choose ri to be the 
distance from the centre to Li . The pattern speed is related to 
ri through the expression 

/50(r) 



Q 



(3) 



where <l)(r) is the potential on the equatorial plane. 

2.2. Galaxy models and free parameters 

In this section we describe the potentials that we use to model 
a barred galaxy. Our basic model will be the one introduced 
by Athanassoula (I1992al l in a thorough study of the orbital 
structure of bars. The axisymmetric component consists of the 
superposition of a disc and a spheroid, whose basic parame- 
ters are determined so that the rotation curve of the galactic 
model has the desired characteristics. The disc is modelled as a 
Kuzmin-Toomre disc (Kuzmin I 19561 Toomre [l963l l of surface 
density 

2X-3/2 



2nr,i 



(4) 



The parameters Vd and set the scales of the velocities and 
radii, respectively. The spheroid is modelled using a density 
distribution of the form 



p{r) = Ph 



-3/2 



(5) 



where pt, and r/, determine its concentration and scale length. 
Spheroids with high concentration have high values of ph and 
small values of rt, the opposite being true for spheroids of low 
concentration. Although we give two separate axisymmetric 
components, it is important to note that, in fact, what matters 
in this study is only the total axisymmetric rotation curve and 
not its decomposition into components. 

The bar component is described by a Ferrers ( 118771 ) ellip- 
soid whose density distribution is described by the expression: 



p^il-m^f m < 1 
m > 1, 



(6) 



where - x^/a^ + y^/b^. The values of a and b determine 
the shape of the bar, a being the length of the semi-major axis, 
which is placed along the x coordinate axis, and b being the 
length of the semi-minor axis. The parameter n measures the 
degree of concentration of the bar. High values of n correspond 
to a high concentration, while a value of n = is the extreme 
case of a constant density bar The parameter po represents the 
bar central density. For these models, the quadrupole moment 
of the bar is given by the expression 

e„, = Mbia^ - b^)/i5 + 2«), 



where Mh is the mass of the bar, equal to 

Mb = 2(2"+^VafeVor(n + l)T(n + 2)/r(2n + 4) 

and r is the gamma function. Throughout this paper we will 
use the following system of units: For the mass unit we take a 
value of IO^Mq, for the length unit a value of 1 kpc and for the 
velocity unit a value of 1 kms"'. Using these values, the unit 
of the Jacobi constant will be 1 km^ s"^. 

For reasons of continuity and to allow in Paper III a com- 
parison of our results with the results of gas flow in these mod- 
els (Athanassoula |1992b] l, we will use the same numerical val- 
ues for the model parameters as in Athanassoula ( |1992al ). We 
will hereafter refer to this model as model A. It has essen- 
tially four free parameters which determine the dynamics in 
the bar region. The axial ratio a/b and the quadrupole moment 
(or mass) of the bar Q„ (or Mb), will determine the strength 
of the bar. The third parameter is the bar angular velocity, or 
pattern speed, determined by the Lagrangian radius ri (Eq.[3]l. 
The last free parameter is the central concentration of the model 
Pc - Pb + Po- As already mentioned, the basic values for the 
free parameters are set as in run 001 of Athanassoula ( |1992a| ): 
a/b = 2.5, rt = 6, Q„, = 4.5 x 10^ pc = 2.4 x lO"*. Then 
a range around each of these values is explored. The axisym- 
metric component is fixed by setting a maximum disc circular 
velocity of 164.204 at r = 20, and r/, is determined by fixing the 
total mass of the spheroid and bar components within r = 10 
to 4.87333 X 10"*, while fixing the combined central density of 
the bar and bulge to p^. The length of the bar is fixed to 5. 
Results on the orbital structure underlying spirals and rings in 
such models will be presented in Sect. 14. II 

So far, we have followed exactly the model of Athanassoula 
( |1992al l. Nevertheless, in this paper we are interested in spi- 
rals and rings which occur beyond or around CR, contrary 
to Athanassoula ( |1992al |1992bl ), who concentrated on orbital 
structure and gas flow within and around the bar region. We 
will, therefore, consider a further option, namely whether the 
rotation curve in the region beyond CR is flat (model F), or 
somewhat rising (model R), or somewhat falling (model D). 
This will be achieved by considering the same axisymmetric 
components as in model A, but with different values of the pa- 
rameters. The different shapes of the rotation curve in the outer 
parts are obtained by giving different maximum disc circular 
velocities at a fixed radius. Thus, model R, with a rising rota- 
tion curve, has a maximum disc circular velocity of 164.204 at 
r - 20. Model F, with a flat rotation curve in the outer parts, 
has a maximum disc circular velocity of 100 at r - 20, while 
model D, with a decreasing rotation curve in the outer parts, 
has a maximum disc circular velocity of 10 at r = 20. The 
spheroid scale length is determined as in model A. In this case, 
we use an inhomogeneous Ferrers bar {n - 1) with axial ratio 
ajb - 2.5 and a semi-major axis a - 5, and a central concen- 
tration Pc - 0.05 X 10"*. The basic values for the remaining free 
parameters are taken to be r/, = 6 and Q„, - 4.5 x 10"*. As 
for the previous model, a range around each of these values is 
explored. Results using these models are presented in Sect. 14.21 

The Ferrers bars are realistic models of bars and have 
been widely used so far in orbital structure studies within and 
in the immediate neighbourhood of bars (e.g. Athanassoula 
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et al. [19831 Pfennig er [T9841 [19871 [T9901 Skokos, Patsis & 
Athanassoula |2002al |2002bl l. They contain parameters with 
physical meaning, such as the bar mass or axial ratio, that can 
be obtained from, or compared to, observations. They have, 
however, one disadvantage, namely that in models using such 
bars the ratio of the non-axisymmetric component of the force 
to the axisymmetric component of the force decreases very 
abruptly beyond a certain radius so that the axisymmetric com- 
ponent dominates in the outer regions. This is of no importance 
if one is interested in the orbital structure or the gas flow in 
the bar region, but in studies like this one, where one is inter- 
ested in the region outside the bar, this may introduce a bias, 
since models with high nonaxisymmetric forces beyond CR 
will not be studied. In order to remedy this, we introduce two 
other models, also often used in the literature, which have an 
ad hoc bar potential, i.e. a potential that is not associated with 
a particular density distribution. Ad hoc models have some dis- 
advantages. They are simple mathematical expressions for the 
potential and do not originate from a realistic density distri- 
bution. This means that the corresponding density distribution 
may have some undesired features, e.g. for very strong nonax- 
isymmetric perturbations the local density could even be neg- 
ative. Furthermore, they do not contain simple parameters that 
can be directly and straightforwardly associated to observable 
quantities, like the bar length, mass, or axial ratio. Most of them 
are of the form eA(r) cos(20), i.e. contain no cos(m6) terms with 
m > 2. This means that the parameter e is associated with the 
mass of the bar and that there is no parameter to regulate its ax- 
ial ratio. Despite all these shortcomings, ad hoc potentials have 
been widely used because they have the important advantage 
of being adaptable to the problem at hand, i.e. with a proper 
choice of the function A{r) one can obtain a potential with the 
desired properties, for example, in our case, a potential with an 
important m - 2 contribution between CR and OLR. 
The first ad hoc potential we use has the form 



1 9 

Q)(r,0) = --ev2^os(26') 



(7)" 



r < a 



r > a. 



(7) 



The parameter ff is a characteristic length scale of the bar po- 
tential and v'o is a characteristic circular velocity. The param- 
eter e is a free parameter related to the bar strength. Dehnen 
( I2OOOI 1 used this potential with n = 3 to study the effect of 
the OLR of the bar of our Galaxy on the local stellar ve- 
locity distribution. Fux (1200 11 1 uses the same bar potential to 
model the Galactic bar and to study the order and chaos in 
the disc. Nevertheless, since the characteristics of the potential 
with n = 3 are very similar to Ferrers potentials, i.e. the force 
decreases very abruptly at large radii, we will use here n = 0.75 
to avoid this. We will couple this bar with an axisymmetric part 
given by Eq.|4]and|5]and we will consider three different slopes 
of the rotation curve in the outer parts, like we did for the pre- 
vious models. In all three models, the length scale of the bar is 
set to ff = 5 and the circular velocity vq is set to 200. The basic 
value for the corotation radius is set to ri = 6.0, and for the bar 
strength to e = 0.15 and a range around each of these values 



is explored. The spheroid mass is fixed at M 1,1,1 = 3x10'*, and 
the scale length is calculated accordingly. In model R', which 
has a rising rotation curve, the disc parameters are fixed so that 
the maximum circular velocity is 164.204 at r = 10. In model 
F', which has a flat rotation curve in the outer parts, we fix the 
values of r,/ and Vd so that the maximum circular velocity is 
164.204 at r = 50. Model D', which has a decreasing rotation 
curve in the outer parts, has a maximum disc circular velocity 
of 164.204 at r = 100. Results using these models are presented 
in Sect. 1431 

We also used in our computations the bar potential given 
by the expression: 



(b{r,d) = eVr(ri - r)cos(26i). 



(8) 



where ri is a characteristic scale length of the bar poten- 
tial, which we will take for the present purposes to be equal 
to 20. The parameter e is related to the bar strength. This 
type of model has already been widely used in studies of 
bar dynamics (e.g. Barbanis & Woltjer 119671 Contopoulos 
& Papayannopoulos 119801 Contopoulos 119811 Athanassoula 
II99OI 1. We will couple this bar with an axisymmetric part given 
by Eq. |4] and |5] We will consider two reference models noted 
with the names of models R" and D", again each with a differ- 
ent slope of the rotation curve in the outer parts. The strength 
parameter is set to e = 100 and the Lagrangian radius is = 6. 
The spheroid mass is Mhui = 3 x 10"*, and the spheroid scale 
length is calculated accordingly. The disc parameters in model 
R" are set so that the maximum disc circular velocity is 164.204 
at r = 30, which corresponds to a rotation curve rising in the 
outer parts, while for model D" they are obtained by fixing a 
maximum circular velocity of 164.204 at r = 10, which corre- 
sponds to a rotation curve falling in the outer parts. 

3. Invariant manifolds and transfer of matter 

In this section, we give a simplified description of the dynam- 
ics around the unstable equilibrium points. Readers interested 
in a more thorough description, in particular of the Lyapunov 
periodic orbits and of the corresponding invariant manifolds, 
are referred to Paper I and to references therein. Here, we also 
define the heteroclinic and homoclinic orbits and describe the 
role they play in the transfer of matter within the galaxy. We 
will use these ideas to compare our results with observations in 
the forthcoming Paper III. 

In Sect. 12. II we mentioned that the equations of motion (Eq. 
|2]l have five equilibrium points, three of which are linearly sta- 
ble, namely L3, L4, and L5, and two unstable, namely Li and 
L2- Around the equilibrium points there exist families of peri- 
odic orbits, e.g. around the central equilibrium point the well- 
known xi family of periodic orbits which is responsible for 
the bar structure. Around each unstable equilibrium point there 
also exists a family of periodic orbits, known as the family of 
Lyapunov orbits (Lvapunov |1949l l, which, at low energy levels, 
are unstable and become stable only at high energies (Skokos et 
al. |2002al l. We are interested in this family only in the range of 
energies where the periodic orbits are unstable. In this range, 
the size of the periodic orbits remains small and they stay in 
the vicinity of the equilibrium point. For a given energy level 
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Fig. 1. Homoclinic (a), heteroclinic (b) and escaping (c) orbits (black thick lines) in the configuration space. In red lines, we 
plot the unstable invariant manifolds associated to the periodic orbits, while in green we plot the corresponding stable invariant 
manifolds. In dashed lines, we give the outline of the bar and, in (b) and (c), we plot the zero velocity curves in dot-dashed Unes. 



within this range, two stable and two unstable sets of asymp- 
totic orbits emanate from the periodic orbit and they are known 
as the stable and unstable invariant manifolds, respectively. We 
denote by the stable invariant manifold associated to the 
periodic orbit y around the equilibrium point L,, / = 1,2. This 
stable invariant manifold is the set of orbits that tends to the pe- 
riodic orbit asymptotically. In the same way we denote by W", 
the unstable invariant manifold associated to the periodic orbit 
y around the equilibrium point L,-, ; - 1,2. This unstable in- 
variant manifold is the set of orbits that departs asymptotically 
from the periodic orbit (i.e. orbits that tend to the Lyapunov or- 
bits when the time tends to minus infinity). Since the invariant 
manifolds extend well beyond the neighbourhood of the equi- 
librium points, they can be responsible for global structures. 

In our planar model, the dynamics take place in a four 
dimensional phase space. In fact, the number of dimensions 
can be reduced by one by fixing the energy to the energy 
level of the Lyapunov periodic orbit. In this three dimensional 
space, the stable and unstable invariant manifolds associated to 
a Lyapunov orbit are sets of asymptotic trajectories that form 
tubes. These tubes are filled and surrounded by a bundle of 
trajectories (see Fig 5 of Paper I, and Gomez et al. 120041 and 
references therein for more details). A good way to visualise 
these tubes and make use of them is by means of Poincare sur- 
faces of section, this is, by drawing the crossings of the tra- 
jectories through a particular plane, or surface in phase space. 
Depending on the purposes of our study, some surfaces will be 
more suitable than others, but the methodology is the same. Let 
us take as an example the surface of section S defined hy y - Q 
with X > Q, this is, we consider the orbits when they cut the 
plane y = having a positive value for the x coordinate. Let 
us consider this surface of section S for the stable and unstable 
invariant manifolds of a Lyapunov orbit around L2 (located in 
the jc < side). Taking initial conditions close to the Lyapunov 
orbit and integrating W", forward in time (resp. W^, backwards 
in time) until the first encounter with S we obtain the simple 
closed curves W";' (resp. W^;') which can be seen in Figs. |2] 
and|3]that will be commented later. In Fig. [T^ we can also see 
the plot of the trajectories until the encounter with S. Although 
the simple closed curves W;;' and w;;' are obtained as the nat- 
ural result of intersecting the manifold tubes with a plane, it 
should be mentioned that further crossings (i.e. w;f and W^f 



with A: > 1) may well not have this simple structure, as can be 
seen in Gidea and Masdemont (I2007I I. 

In the selected example, W";' and Wy'} are represented in 
(x, x) coordinates. It is important to note that a pair (x, i) in S 
defines an orbit in a unique way, since y - and y is obtained 
from the energy level under study. By the definition of an in- 
variant manifold, a point in W"'^^ n Wy^ represents a trajectory 
asymptotic to the Lyapunov orbit y around L2 both forward 
and backward in time. It is called an homoclinic orbit. In gen- 
eral, homoclinic orbits correspond to asymptotic trajectories, 
iff, such that iff € Wy. n Wy,, i - 1,2. Thus, a homoclinic or- 
bit departs asymptotically from the unstable Lyapunov periodic 
orbit y around L, and returns asymptotically to it (see Fig.[T^). 
Heteroclinic orbits, on the other hand, are defined as asymp- 
totic trajectories, iff', such that tf/' e W" n Wy., i + j, i, j — 1,2. 
Thus, a heteroclinic orbit departs asymptotically from the peri- 
odic orbit y around L, and approaches asymptotically the corre- 
sponding Lyapunov periodic orbit with the same energy around 
the Lagrangian point at the opposite end of the bar Lj, i + j (see 
Fig. [TJ); a suitable surface of section for this computation can 
be the plane x - 0). 

In our computations we also consider the overlap area cor- 
responding to homoclinic orbits as the area resulting from the 
intersection of the interior regions of W"'' and Wy^ , i - 1,2 
(see Fig. |3]l. Analogously, we define the overlap area corre- 
sponding to heteroclinic orbits as the area resulting from the in- 
tersection of the interior regions of W"'' and W^' ', / + j, i, j — 
1,2. Wy.^ and Wyf for a given model and energy may or may 
not intersect and, as the parameters of the model change, they 
can approach each other or move away. For instance. Fig. |2] 
shows the curves W;'' and Wyl^ on the plane S for our model 
F, introduced in Sect. 14.21 for four different values of the 
quadrupole moment Q,„. We note that models with nonzero in- 
tersection of Wy'.^ and Wy^ are not associated to isolated values 
of the given free parameter (in our example Qm), but they de- 
fine a range. So we will find trajectories that depart asymptoti- 
cally from the periodic orbit, follow W", and do not intersect the 
corresponding stable invariant manifold in phase space. These 
trajectories initially spiral out from the region of the unstable 
periodic orbit (see Fig. [1]:) and we refer to them as escaping 
trajectories. 
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We will argue in Paper III that these three types of orbits 
- namely, the homoclinic, the heteroclinic, and the escaping 
orbits - become the backbone of ringed structures and of spi- 
ral arms observed in disc galaxies and we will follow how the 
overlap area gives a measure of these behaviours. 

4. Results 

The next subsections are devoted to the results obtained when 
we vary the parameters of the models introduced in Sect. 12.21 In 
order to best see the influence of each parameter separately, we 
make families of models in which only one of the free parame- 
ters is varied, while the others are kept fixed. We start by study- 
ing the effects of the main four parameters of model A. We then 
study the effect of having reference models with a Ferrers bar 
potential and rising, flat, or falling rotation curves in the outer 
parts. Finally, we study the effects of having bar potentials with 
a larger influence in the outer parts. 

4.1. Model A. The effect of the main four parameters 

In this section we study the effect of the variation of the pa- 
rameters of model A on the shape of the invariant manifolds. 
We consider both homogeneous (n = 0) and inhomogeneous 
(n = 1) Ferrers bar models (Eq. |6]l. The main four parame- 
ters are the axial ratio of the bar, a/b, the Lagrangian radius, 
ri, the quadrupole moment, Q,,,, and the central concentration, 
Pc- For all four, we will take the same range of values as in 
Athanassoula (|1992al l, since it is wide enough to make sure we 
cover all relevant values. According to Kormendy ([1982]), typ- 
ical bar axial ratios range from 2.5:1 to 5:1, while ovals have 
considerably smaller axial ratios. We thus choose values of the 
axial ratio within the range a/b - 1.2 - 6. As the rings and 
spirals generally emanate from the ends of bars, the unstable 
Lagrangian points should be placed in the vicinity of the bar 
end points, which are at x = ±5. We let the parameter vary 
accordingly within the range - 3.5 - 8.5, i.e. we consider 
Lagrangian points placed from well within the bar up to well 
beyond the bar end points. The values for the quadrupole mo- 
ment are chosen within the range Q„, - 0.1 x 10"* - 12 x 10^, 
covering all the range from weak bars or ovals to very strong 
bars. Finally, the central concentration is related to the presence 
of the spheroid. In this set of models the central concentration 
is varied within the range pc = 0.02 x 10"* - 3.6 x 10"*. For a 
given value of the bar central density, higher values of pc in- 
dicate denser bulges. We used two values of the concentration 
index of the bar, namely n = and n - \. For each of these 
indexes, we make families of models in which just one of the 
free parameters is varied, while the others are kept fixed. We 
use these families to study the influence of each of the free pa- 
rameters on the global shape of the invariant manifolds. 

In Fig.|5]we show the effect of the four free parameters on 
the shape of the invariant manifolds for model A with n - 0. 
In each panel, we plot the outline of the bar and, for a given 
energy level, the unstable invariant manifolds associated to the 
Lyapunov orbits of both Li and L2 equilibrium points, inte- 
grated until they perform half a revolution. For initial condi- 
tions of the invariant manifolds very near the Lyapunov orbit. 




n=l _ 
n=0 - - 
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Fig. 4. Rotation curves corresponding to model A for « = 
(dashed line) and n = 1 (solid line). 

this time corresponds in all cases to approximately 3 bar rota- 
tions, i.e. to ~ 0.25 Gyr In columns 1 to 4 we show the effect of 
the bar axial ratio (first column), of the Lagrangian radius (sec- 
ond column), of the quadrupole moment (third column) and of 
the central concentration (fourth column). In all columns the 
numerical values of the parameters increase from top to bot- 
tom. Results for n - \ are very similar, so we do not show 
them here. 

Fig. |5] shows that the central concentration has hardly any 
influence on the shape of the invariant manifolds, which corre- 
sponds in all cases to the morphology of an rRi ringed galaxy. 
From the ends of the bar emanate branches that curl around the 
bar until the opposite end, forming two rings. One ring is close 
to the bar and is elongated along it, but not far from circular, 
i.e. has the location and properties of observed inner rings. The 
other ring is at larger radii and is elongated perpendicularly to 
the bar major axis, i.e. it has the location and properties of an 
outer ring Ri. This description also holds for all values of the 
bar axial ratio. Varying this quantity has hardly any influence 
on the outline of the outer branches of the invariant manifolds, 
but does have some, albeit small, on the inner ones. Namely, 
for ajb - 1.2, the inner ring is near-circular, while for larger 
values it is diamond shaped. 

On the contrary, varying the value of introduces different 
morphologies. When the Lagrangian points are well within the 
bar, the inner ring is also within the bar, which is not realistic. 
The outer ring is not closed, i.e. we have a spiral. However, this 
spiral also is not realistic, because it does not emanate from 
the ends of the bar For values of the Lagrangian radius beyond 
ri - 5, the outer branches of the invariant manifolds tend to 
the opposite ends of the bar, and when they close they form an 
outer ring. Thus, the morphology is again that of an rR\ ringed 
galaxy. As ri increases, the inner ring increases in size and 
becomes more circular. Already for ri -6.5 the inner ring size 
is considerably larger than the bar, contrary to what is observed. 
For the highest of the considered values, the inner ring forms a 
nearly circular structure with about double the size of the bar 

Finally, the variation of the bar strength also gives different 
shapes of the ringed structures. In the case of weak bars, i.e. 
small Qm values, the inner ring is nearly circular with a diam- 
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Fig. 2. Evolution of the curves W";' (black solid lines) and Wy} (red dashed line) on the phase space xx for a particular galactic 
model and four values of 2m (given in the upper left corner of each panel in units of 10^). 




Fig. 3. Definition of the overlap area corresponding to homoclinic orbits. The curves W";^ (in red solid hne) and W^;' (in green 
dashed line) intersect on the plane (x, x). The four black points correspond to homocUnic orbits, while the hatched area is the 
overlap area corresponding to homoclinic orbits. 



eter about the size of the bar However, as will be discussed 
in paper III, the outline is not that of an rR\ ringed galaxy. As 
the quadrupole moment increases and the bar gets stronger, the 
inner ring becomes more elongated along the bar until, for the 
case of strong bars, it is deformed in shape. The outer ring gets 
bigger in size with increasing 2,„, so that the ratio between the 
ring diameters increases. In the case of strong bars, both the 
inner and outer rings become asymmetric with respect to the 



major and minor axes of the bar and their major and minor di- 
ameters do not coincide any more with the x and y axes. 

We can therefore conclude that the shape of the invariant 
manifolds is not sensitive to the variation of the axial ratio 
and of the central density. On the other hand, when we vary 
the Lagrangian radius or the quadrupole moment, we obtain 
richer structures. These differences are related to the values of 
the effective potential of the models. This can be seen in Fig.|6] 
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Fig. 5. Ring and spiral structures in four sequences of models with n - and various values of the four main parameters (given 
in the upper corners of each panel). For the remaining parameters, see text. First column: effect of the variation of the axial ratio 
of the bar, a/b; Second column: effect of the variation of the Lagrangian radius, ri; Third column: effect of the variation of the 
quadrupole moment, Q,„ (in units of 10"*); Fourth column: effect of the variation of the central density, pc (in units of 10*). In all 
panels, we plot on the x - y plane the outline of the bar (thin solid line) and the unstable invariant manifolds associated to the 
Lyapunov orbits for a given energy level. 



where we plot the effective potential of the models along the bar 
major axis. First, we note that for radii smaller than the length 
of the bar, the effective potential shows some differences in all 
the cases. This means that all four free parameters will have 
an effect in the interior structure of the bar and is in agreement 
with the fact that both the orbital structure and the gas flow in 
the bar region vary considerably as the values of the four main 
parameters are changed (Athanassoula |1992al |1992bl l. Here, 



however, we are interested in the spirals and rings which occur 
in a region from the ends of the bar outwards. For radii larger 
than the length of the bar, the effective potential depends little 
on the value of the axial ratio, or of the central concentration. 
This is illustrated in Fig. and Fig. |6}; for homogeneous and 
inhomogeneous bars, respectively. On the other hand, when we 
vary the Lagrangian radius or the quadrupole moment, the ef- 
fective potential is considerably influenced. This is shown in 
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Fig. 7. Rotation curves corresponding to models R (solid line), 
F (dashed line) and D (dot-dashed line). Note that all models 
have the same rotation curve slope in the inner parts. 

Fig. |6j3 and Fig. |6]i, again for homogeneous and inhomoge- 
neous bars, respectively. This explains why only the variation 
of the Lagrangian radius and of the quadrupole moment has an 
effect on the global shape of the invariant manifolds. 

4.2. Models R, F and D. The effect of having a rising, 
flat, or falling rotation curve with a Ferrers ellipsoid 

In general, rotation curves in the outer parts of disc galaxies 
are flat (Bosma[198T). However, many cases of slightly rising 
or slightly falling rotation curves are known. We thus want to 
study the possible influences of these different rotation curves 
shapes on the morphological structure of our models. 

In this section we use reference models R, F and D, each 
having a different slope of the rotation curve in the outer parts. 
The rotation curves of these models are shown in Fig.|7] Note 
that in all cases the slope of the inner parts is the same, while in 
the outer parts it is rising (model R), flat (model F), or decreas- 
ing (model D). For each of these reference models we generate 
families of models, where we vary only one of the free param- 
eters. Since in Sect. 14. II we showed that the whole range of two 
of the main parameters corresponds to the same morphology 
and dynamics beyond the bar region, and explained the reason 
for this, here we will study only the effect of the variation of 
the two crucial parameters, namely the values of ri and Q„. 
The range of variation of these parameters will be the same as 
in Sect.gJl i.e. = 3.5 - 8.5 and Q„r = 0.5 x lO"* - 12 x 10^ 
while the values of pc and a/b are kept fixed to 0.05 x 10"* and 
2.5, respectively. 

In Fig. [8] we show the effect of the variation of the 
Lagrangian radius on the shape of the invariant manifolds. In all 
columns we increase the value of ri from top to bottom, while 
the value of the quadrupole moment is fixed to Q„, = 4.5 x 10"*. 
We note that, for rotation curves which in the outer parts are 
rising or flat, the shape of the invariant manifolds presents the 
same behaviour as in model A (see previous subsection). Thus, 
for small values we have spiral arms, while for large val- 
ues we have outer rings. An inner ring is present in all cases 
and its size increases with r^. It is smaller than the bar, i.e. un- 
realistic, for ri < a and much larger than the bar, i.e. again 



unrealistic, for very large ri values, while realistic sizes are 
only found for an intermediate range of values. In the case of 
falling rotation curves, the evolution is similar. However, when 
the Lagrangian radius is close to the bar length, the invariant 
manifolds cross each other, forming an R2 ring. 

Along each row we observe the effect of changing the slope 
of the rotation curve at the outer radii, while keeping constant 
the value of all bar parameters, including the corotation radius. 
In the upper row, where the Lagrangian points are inside the 
bar, the spiral arms open as the slope of the rotation curve de- 
creases. A case where corotation is near the ends of the bar is 
shown in the second row. A rising rotation curve gives an rR\ 
ring, while a flat rotation curve gives a Ri pseudo-ring. In the 
specific case of a falling rotation curve, the outer branches of 
the invariant manifolds are long and open, finally crossing each 
other, thus forming an R2 ring. For values of the corotation ra- 
dius beyond the end of the bar, the shape corresponds to that of 
an rR\ ring galaxy. 

In Fig. |9l we show the effect of the variation of the 
quadrupole moment on the shape of the invariant manifolds. 
In all columns we increase the value of Q,„ from top to bottom, 
while the value of the Lagrangian radius is fixed to = 6. In 
the case of reference models with rising or flat rotation curves 
(first and second columns), the behaviour of the invariant man- 
ifolds is the same as in model A. For a falling rotation curve 
and a weak bar, we still have an rR\ ring structure. As we in- 
crease the bar strength, first the inner ring becomes distorted, 
then the major axis of the outer ring ceases to be perpendicular 
to the bar major axis, and finally, with strong bars, the outer 
invariant manifolds either do not close and form spiral arms, 
or they are long and cross each other, thus forming R2 rings. 
As in the previous figure, along each row we observe the effect 
of changing the slope of the rotation curve in the outer parts, 
while keeping constant the bar quadrupole moment. For weak 
and intermediate strength bars, the shape of the rotation curve 
does not influence the morphology. For strong bars, however, 
we obtain different structures. A rising rotation curve gives an 
rR\ ring, while a flat rotation curve gives an asymmetric ring. 
In the case of falling rotation curves, the outer branches of the 
invariant manifolds are open and form spiral arms and Ro rings. 

4.3. Models R', F' and D'. Other bar models 

As previously mentioned, Ferrers ellipsoids do not have a great 
influence in the outer region, which is where the outer branches 
of the invariant manifolds are found. Since we argue in this pa- 
per that invariant manifolds are responsible for the rings and 
spiral structures, we considered also bar potentials of either 
of Dehnen type (Eq. |7]i, or of Barbanis-Woltjer type (Eq. [8]i. 
The results of the two models are, however, very similar, so for 
the sake of brevity we discuss only the results for the Dehnen 
type bar here. In this case we use models R', F' and D', again 
each with a different slope of the rotation curve in the outer 
parts. The rotation curves of these models are shown in Fig. 
[TOl As in the previous sections, we make families of models, 
varying one of the parameters, while keeping the rest constant. 
Since these models are ad hoc and they do not correspond to 
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Fig. 6. Effective potential along the bar semi-major axis for model A with n = and n - 1 . (a) Models with homogeneous (n - 0) 
bars with different values of the bar axial ratio, ajb, and of the central density pc (in units of 10'*). (b) Models with homogeneous 
bars (« = 0) with different values of the Lagrangian radius, r^, and of the quadrupole moment, Q„, (in units of 10"*). (c) Same 
as in (a) but with inhomogeneous bars (n = 1). (d) Same as in (b) but with inhomogeneous bars (n — 1). The solid vertical line 
marks the end of the bar 
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Fig. 10. Rotation curves corresponding to models R' (solid 
line), F' (dashed line) and D' (dot-dashed line). Note that all 
models have the same rotation curve slope in the inner parts. 

any simple and realistic density, choosing the parameters is less 
easy, since there are less constraints from observations. Based 
on the results of Sect. 14.11 we use as free parameters only the 
Lagrangian radius, ri, and e, which determines the bar strength. 
The values of ri are taken within the range ri - 3.5 - 8.5, as in 
the previous sections. The values of e are taken within the range 
e - 0.01 - 0.3, which corresponds a range of 5 to 26% for the 
ratio of the non-axisymmetric component to the axisymmet- 
ric component of the force in the x-direction and at the radius 
r = 16. 



In Fig.[TT]we show the effect of the Lagrangian radius on 
the shape of the invariant manifolds. In all columns we increase 
the value of ri from top to bottom, while the strength param- 
eter € is fixed to 0.15. For model R', i.e. for a rotation curve 
which is slightly increasing, we start at the smallest ri with 
a spiral. As the value of is increased the spiral gets more 
tightly wound and then closes forming a somewhat asymmet- 
ric outer ring. These changes are accompanied by an increase 
of the size of the inner ring, which, nevertheless, is in all cases 
quite asymmetric. For model D', i.e. for a somewhat decreasing 
rotation curve, at the smallest ri value we have a rather open 
spiral. As the value of is increased, the extent of the spi- 
ral increases considerably and its pitch angle increases slightly. 
For the largest ri value, - 8.5, the outer branches of the 
stable and unstable invariant manifolds delineate well the loci 
of both types of outer rings, namely Ri and R2, thus forming 
an R1R2 ring. The inner ring also increases in size as the value 
of ri increases, while staying of roughly the same shape. The 
position angle leads that of the bar by an angle which, as we in- 
crease kl, increases approximately from 20 to 36 degrees. The 
sequence with the flat rotation curve, F', is the richest in mor- 
phological types. It starts with a rather open spiral, as the D' 
sequence, but, contrary to the D', the shape of the invariant 
manifold branches changes considerably. Thus, for the largest 
ri value they form asymmetric outer and inner rings, similar in 
outline to those formed for model R' and - 4.5. In between 
these two, the invariant manifolds give us two different mor- 
phologies. For ri - 4.5, the outer branches of the stable and 
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Fig. 8. Rings and spirals for three sequences of models with reference models R (left), F (middle) and D (right) where we vary 
the Lagrangian radius, r£. This value increases from top to bottom and is given in the upper right corner of each panel. The value 
of the quadrupole moment is fixed to Qm = 4.5 x 10"*. In all panels, we plot on the x - y plane the outline of the bar (thin solid 
line), the two unstable Lyapunov orbits associated to the unstable equilibrium points L\ and L2 (white solid line), and the unstable 
invariant manifolds associated to the Lyapunov orbits, all for a given energy level. 



unstable invariant manifolds outhne the shape of R1R2 rings, as 
for model D' and = 8.5, while the inner ring is deformed in 
shape. For = 5.5, the outer branches of the unstable invariant 
manifolds evolve until they cross each other forming an outer 
ring whose principal axis is parallel to the bar major axis, thus 
forming an R2 outer ring. 



In each row we observe the efi'ect of keeping constant the 
value of corotation radius and of all other bar parameters, while 
changing the slope of the rotation curve at outer radii. When 
the Lagrangian radius is smaller than the bar length, the spiral 
arms open as the slope of the rotation curve decreases. The case 
where corotation is approximately the value of the bar length 
scale is shown in the second row. A rising rotation curve gives 
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Fig. 9. Same as in Fig. [8] for three sequences of models with reference models R (left), F (middle) and D (right) and different 
values of the quadrupole moment, Q,„, given in units of 10"* in the upper right corner of each panel. The value of the Lagrangian 
radius is fixed to ri = 6. 



a broken and misaligned rRi ring, while a flat rotation curve 
gives RiRi rings, and a falling rotation curve gives spiral arms. 
The tendency is similar if the values of the Lagrangian radius 
are somewhat larger than the bar length scale. For a rising ro- 
tation curve, the structure of the galaxy is that of an rRi ring. 
However, if the rotation curve is flat, we obtain R2 rings, while 
if the rotation curve is falling, we obtain spiral arms. As we in- 
crease the value of the Lagrangian radius further, the invariant 
manifolds tend to approach the opposite ends finally forming 



outer rings misaligned with respect to the bar. The inner rings, 
in turn, tend to open as the Lagrangian radius increases. In the 
case of a falling rotation curve, the invariant manifolds form 
spiral arms. 

In Fig. [12] we show the effect of the bar strength on the 
shape of the invariant manifolds. In all columns we increase the 
value of £ from top to bottom, while the value of the Lagrangian 
radius is fixed to = 6. In the case of reference models with 
a rising rotation curve, the invariant manifolds tend to open. 
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Thus, we have a behaviour opposite to the one we found when 
we increased the value of the Lagrangian radius. That is, for 
low values of the parameter e, the galaxy morphology is that 
of an rRi ring, while as the bar strength increases, the branches 
of the invariant manifolds tend to open and form spiral arms. 
The inner ring is, for low values of e, symmetric. As the bar 
strength increases, it becomes more eccentric and finally asym- 
metric. For a flat rotation curve and a weak bar, we initially 
have an rRi ring structure. As we increase the bar strength, the 
inner ring becomes distorted, while the major axis of the outer 
ring ceases to be perpendicular to the bar major axis. Increasing 
the bar strength, we find spirals. A further increase of the bar 
strength brings a morphology of the R1R2 type (see middle 
panel for e = 0.2) and then (see middle panel for e - 0.25) 
a morphology of the R2 type. Finally, with strong bars the outer 
invariant manifolds do not close, thus forming spiral arms. For 
a falling rotation curve, the behaviour is similar to the F' case 
but the invariant manifolds become open for weaker bars. As 
the bar strength increases, the outer branches of the invariant 
manifolds open forming R1R2 rings and spiral arms, but the in- 
ner branches will appear distorted. 

In each row, we observe the effect of changing the slope of 
the rotation curve in the outer parts. For weak bars, the invari- 
ant manifolds tend to open as the slope of the rotation curve 
decreases. Note that this behaviour is similar to the one ob- 
tained for keeping constant the value of corotation radius. We 
obtain R1R2 rings for weak bars and falling rotation curves or 
intermediate strength bars and flat rotation curves. In the case 
of strong bars, the morphologies obtained are spiral arms and 
R2 rings. 

5. Summary and discussion 

In this paper we developed an idea originally proposed in Paper 
I, namely that spirals and rings in barred galaxies can be ex- 
plained in a common dynamical framework. This is based on 
the dynamics of the unstable equilibrium points in the rotating 
bar potential and, more specifically, on the invariant manifolds 
associated to the unstable Lyapunov periodic orbits around 
these equilibrium points. These unstable Lyapunov orbits are 
simple orbits of elliptical-like shape, encircling one of the un- 
stable equilibrium points and staying always in its immediate 
vicinity. The dynamics of the corresponding manifolds have 
been well studied by orbital structure theory and we have drawn 
heavily from that field for our purposes. It is not, however, nec- 
essary to master the intricacies of the definition and dynamics 
of invariant manifolds in order to follow our work. The reader 
may simply think of a set (a bundle) of orbits emanating from 
the vicinity of the Lyapunov orbits, i.e. from nearby positions 
and velocities. 

We investigated different barred galaxy models, to make 
sure that none of our principal results are model dependent 
and also in order to give ourselves a broad base to allow us, 
in Paper 111, to discuss globally the application of our results 
to real galaxies. Model A is identical to that introduced by 
Athanassoula (1992a). Models R, F and D are similar, but al- 
low for rotation curves which can be slightly rising (R), flat 
(F), or decreasing (D). Finally, Models R', F' and D' (or R" 



and D") are similar to R, F and D, but have a Dehnen-type (or 
a Barbanis-Woltjer type) bar 

Our main result is that the loci outlined by the invariant 
manifolds and by the orbits associated with them can repro- 
duce, for appropriate values of the model parameters, all the 
morphologies observed in real galaxies. Thus we can obtain 
two spiral arms emanating from the ends of the bar, or rings. 
We obtain inner rings as well as all the varieties of outer rings, 
namely Ri, R2 and RiR2- This is shown in Fig. [13] where we 
show one typical example of each kind of morphology. These 
include an rRi morphology (Fig. [T3h) as in NGC 2665 (Buta 
& Crocker tl991jl, an rR2 one (Fig. as in ESO 325 - 28 
(Buta & Crocker[T99B, mRiR2 one (Fig.[T3};) as in NGC 3081 
(Buta & Purcell [T998l Buta, Byrd & Freeman l2004l l. and a spi- 
ral (Fig.injl) as in NGC 1365 (Jorsater & Moorsel [T995] l. It is 
important to stress that all these types were obtained with the 
same dynamics. Indeed, all that differs between these four ex- 
amples is the values of the model parameters and the bar model 
potential. 

Model A has four basic free parameters describing the bar 
dynamics. We show that two of these, the bar axial ratio and the 
central concentration of the model, hardly influence the loci 
of the invariant manifolds, while the two others, namely the 
Lagrangian radius and the bar quadrupole mass, have a con- 
siderable effect. We explain this by looking at the influence of 
these four parameters on the effective potential of the system. 
The main type of morphology found in model A is the rRi type. 
Allowing for different slopes of the rotation curve and/or dif- 
ferent bar models, we introduce the remaining morphological 
varieties. 

In this paper we considered a very large number of models 
and of model parameters and calculated in all cases the appro- 
priate invariant manifolds in order to find the resultant mor- 
phology. This can also be applied to real galaxies. Using an 
image of the galaxy in the near infrared and with an estimate of 
the thickness and of the mass-to-light ratio, one can obtain the 
bar potential. Including information from the rotation curve, 
one can obtain the total potential in which the manifolds and 
orbits can be calculated. A similar procedure has been already 
applied e.g. by Lindblad, Lindblad & Athanassoula ( 1996l l to 
study the gas flow in NGC 1365 and by Byrd, Buta & Freeman 
(.2006J to study the rings of NGC 3081. The latter work in 
particular studies the shape of the inner and outer rings both 
analytically and with the help of simulations. Of course, as is 
the case for all modelling of observations, e.g. concerning the 
mass-to-light ratios and the thickness, this procedure relies on a 
few approximations and assumptions. Yet it has revealed a lot 
of important information on bars and spirals (e.g. Sanders & 
Tubbs lT980l Duval & Athanassoula [Mil Lindblad, Lindblad 
& Athanassoula 119961 Patsis, Athanassoula & QuiUen 119971 
Weiner, Sellwood & Williams |200r, Kranz, Slyz & Rix [2 0031 
Perez, Fux & Freeman l2004l Byrd, Buta & Freeman l2006l l and 
would thus be well worth pursuing for the model presented 
here. 

This paper leaves two major points unanswered. The first 
one is why a given set of model parameters gives a given mor- 
phology? Our work so far shows that this is not a random pro- 
cess, but does not explain the physics behind it. Thus we are not 
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Fig. 11. Rings and spirals for three sequences of models with reference models R' (left), F' (middle) and D' (right), where we 
vary the Lagrangian radius, r/.. This value increases from top to bottom and is given in the upper right corner of each panel. The 
value of the strength parameter is fixed to e = 0.15. In all panels, we plot on the x - y plane the two unstable Lyapunov orbits 
associated to the unstable equilibrium points Li and L2, and the unstable invariant manifolds associated to the Lyapunov orbits, 
all for a given energy level. For model F' and = 4.5 and model D' and = 8.5, we plot the corresponding stable invariant 
manifolds. 



able, at this point, to foresee what morphology will result from 
each particular bar model and set of specific parameters with- 
out calculating the necessary invariant manifolds. The second 
unanswered point is to what extent the results found here are 
relevant to real galaxies. We have of course shown that the loci 
of the appropriate invariant manifolds can reproduce all the ob- 
served morphologies. This is very encouraging, but does not. 



on its own, assure us that we have found here the theory that 
explains the formation of ring and spiral structures in barred 
galaxies. Several more points need to be considered before this 
conclusion can be reached. These yet unanswered two points 
will be addressed in Paper III. 
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Fig. 12. Same as in Fig.[TT]for three sequences of models with reference models R' (left), F' (middle) and D' (right) and different 
values of the bar strength parameter, e. The value of the Lagrangian radius is fixed to - 6. For model F' and e - 0.2 and model 
D' and e = 0.1, we plot the corresponding stable invariant manifolds. 
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